Phonon-Induced Backscattering in Helical Edge States 
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A single pair of helical edge states as realized at the boundary of a quantum spin Hall insulator is 
known to be robust against elastic single particle backscattering as long as time reversal symmetry 
is preserved. However, there is no symmetry preventing inelastic backscattering as brought about 
by phonons in the presence of Rashba spin orbit coupling. In this letter, we show that the quantized 
conductivity of a single channel of helical Dirac electrons is protected even against this inelastic 
mechanism to leading order. We further demonstrate that this result remains valid when Coulomb 
interaction is included in the framework of a helical Tomonaga Luttinger liquid. 

PACS numbers: 71.10.Pm,72.15.Nj,85.75.-d 



During recent years, great interest has been attracted 
by the theoretical prediction p] [2] and experimental dis- 
covery [3] of the quantum spin Hall (QSH) effect. The 
QSH phase is a two dimensional realization of a topo- 
logical insulator (TI), a class of time reversal symme- 
try (TRS) preserving phases that differ essentially from 
trivial atomic insulators by a Z2 topological invariant 
[3HS]- Besides the enormous conceptual depth of topo- 
logical states of matter, TI phases are also considered 
promising candidates as to future applications in nano- 
electronics. This is due to their topologically protected 
transport properties which might be exploited for high 
precision spintronics devices operating at low power con- 
sumption. Therefore, modelling the QSH effect under 
experimentally relevant conditions is crucial to test the 
practical limitations of these protected features. 



As far as the robustness of the topological protection 
is concerned the QSH effect is fundamentally different 
from the integer quantum Hall (IQH) effect. For the TRS 
breaking (IQH) insulator [THH] topological order leads to 
a quantization of conductivity to impressive accuracy. In 
the IQH regime, edge state transport is chiral, mean- 
ing that the density of states for subgap backscattering 
vanishes which excludes such processes by simple phase 
space arguments. In contrast, in the TRS preserving 
QSH phase a single pair of helical edge states induced 
by bulk boundary correspondence is supported at the 
edge of the QSH bar. This means that both right- and 
lcftmovers exist at a single edge. However, states of op- 
posite direction of motion are Kramers partners due to 
TRS. The well known topological protection of a single 
pair of helical edge states against backscattering in this 
scenario can be mathematically illustrated by the follow- 
ing simple argument [10]. Let |</>) and — T\<j>) be 
Kramers partners. Then as long as H is a TRS preserv- 



ing Hamiltonian 

(1>\H\<j>) ={4>WW = (T<p\TH\iP) = 

(iP\HT\tP) = Mht 2 \4>) = (1) 

i.e. the matrix element for scattering between the 
Kramers partners vanishes. Note that argument ([I]) of 
protection relies on two fundamental constraints: First, 
only single electron processes are considered. Second, 
since Kramers partners are degenerate states, it only 
precludes elastic backscattering. Within the validity of 
these restrictions extensive studies of the helical Tomon- 
aga Luttinger liquid (hTLL) jTQl [TT] representing a sin- 
gle pair of helical edge states have shown that Anderson 
localization is avoided [1] in the presence of TRS pre- 
serving disorder and that TRS breaking magnetic impu- 
rities can open a gap in these systems [12 . Furthermore, 
interedge backscattering can occur if the QSH sample 
is locally narrowed down to a quantum point contact 
[T3HT7] or if two QSH bars are brought close to each 
other [18]. In general, backscattering at a single heli- 
cal edge requires spin flip processes. In realistic setups, 
these are induced by Rashba spin orbit coupling (SOC) 
originating from unavoidable potential fluctuations. Pre- 
serving TRS, Rashba SOC cannot cause single electron 
elastic backscattering, though. However, relaxing the 
single electron processes constraint by additionally in- 
cluding Coulomb interaction, two electron backscattering 
processes have been shown to arise in these systems [TH] ■ 
Such backscattering terms are well known to be allowed 
by TRS [rO- 
Under realistic experimental conditions finite tempera- 
ture and bias voltage also imply the presence of phonons, 
i.e. inelastic processes that undermine the second con- 
straint for the validity of argument 0. It is thus of 
crucial importance to investigate the influence of this dis- 
sipative mechanism on the topological protection. Here, 
we show two important results of the helical edge states in 
presence of two TRS preserving perturbations: Rashba 
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FIG. 1. Pair of helical edge states with two terminals and a 
Rashba impurity coupling the two channels. The wavy line 
illustrates the presence of electron phonon coupling in the 
system. 

SOC and electron phonon coupling. First, we demon- 
strate that in this scenario there is no strict protection 
against inelastic single electron backscattering. Second, 
we find that for helical Dirac fermions the leading or- 
der contribution of this mechanism vanishes, support- 
ing the protection for practical purposes. We further 
demonstrate how this additional robustness fully sur- 
vives in presence of Coulomb interaction, i.e. in a hTLL 
with electron phonon coupling and Rashba SOC. In a 
nonequilibrium transport calculation for the hTLL, we 
take the electron phonon coupling into account exactly 
by integrating out the phonons using a Keldysh con- 
tour path integral representation of the generating func- 
tional. Our analysis is relevant for any realization of 
the hTLL as a one dimensional system. Model with- 
out Coulomb interaction — We investigate a single pair 
of helical Dirac fermions coupled linearly to longitudi- 
nal acoustic phonons. The two species of electrons are 
coupled via Rashba SOC (see Fig. [I]). In most parts of 
this work we will have a sharp impurity-like scattering 
potential in mind which brings about momentum trans- 
fer on the order of 2kp, where kp is the Fermi wave 
vector. We represent the helical fermionic fields as a 
spinor * = (* m ,^> L ;) T = (* + , <F_) T . The free elec- 
tron Hamiltonian then reads 

H u = J dx&(x)pa z $(x) 

where p = —id x is the momentum operator an o~ z is 
a Pauli matrix in spin space. The two most relevant 
Rashba-terms induced by a spatially dependent electric 
field in z-direction are given by [HI [20] 

Hr= \I dx ^^ (M X )>P} + {a 3 {x),p 3 }) <J y ^{x) 

Note that only odd powers of p are allowed by TRS. Elec- 
tron phonon coupling to the displacement field d of longi- 
tudinal acoustic phonons is modelled by the Hamiltonian 

EP 

H cp = X J dx^ (x)a Q ^(x)d x d(x) 

with all dimensionful constants absorbed into A. The free 
phonon dynamics is governed by 

B p = X -ldx [(Tl d (x)) 2 +c 2 (d x d(x)) 2 ] 



where c is the acoustic phonon velocity in units of the 
electronic Fermi velocity and IT^ is the conjugate mo- 
mentum of d. We model the phonons for a strictly ID 
system which corresponds to an in transverse direction 
perfectly localized edge state. Later on, we will see that 
our key results do not critically depend on the details of 
the phonon model. The total Hamiltonian of our setup 
is then given by 

H = -ffhi + H p + Hr + H cp = H + Hi (2) 

where Hq = ify + Hp is the free Hamiltonian whereas 
Hi = Hu + H cp encompasses the coupling terms. 

Inelastic backscattering — We will now demonstrate 
how the combination of Hu and H cp will in principle 
be able to cause single electron backscattering at a sin- 
gle edge of a QSH insulator. As observed above, since 
Hu is TRS preserving it cannot cause elastic single elec- 
tron backscattering. Due to its offdiagonal structure in 
spin space it couples opposite spins though. In contrast, 
H cp does not mix different spin species but can bring 
about energy dissipation by virtue of energy transfer from 
the electronic degrees of freedom to phonons. Thus the 
second order in Hi diagrams shown in Figure [2] which 
are first order in H ep and in Hu cause backscattering at 
finite bias (see Fig. [3| if their contribution does not van- 
ish for momentum transfer pi — pf s» 2k p. We consider 




FIG. 2. Tree diagrams for lowest order backscattering. 
Dashed line with square represents the external Rashba po- 
tential. Wavy line denotes the phonon propagator. 

scattering between a right mover \pf) and a left mover 
\pj,q p h) with an additional phonon. Up to second order 




FIG. 3. Voltage configuration to pass a spin up current from 
the left to the right. The bias V opens an energy window for 
inelastic phonon scattering. 
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the corresponding scattering matrix element Ma can be 
written as 

M i{ =( P j,q ph \H I G H I \ P +) 

where Go is the free propagator corresponding to Hq. In- 
terestingly, the lowest order contribution of the Rashba 
term linear in momentum associated with ot\ vanishes 
due to a nontrivial destructive quantum interference of 
the two contributing diagrams which stems from the lin- 
earity of both H^i and the ai-Rashba term. To show that 
this is not due to any fundamental symmetry like TRS wc 
calculate the same matrix element for the Rashba term 
associated with 123 which yields 



|M if | 2 = 



A 2 c 

1 

16tt 



*3 [Qph+Pf -Pi 



\Qph\ 



where 0,3 (k) is the Fourier transform of 0:3 (x). An analyt- 
ical Fermi's golden rule calculation at zero temperature 
for a (5-shaped Rashba impurity 0:3 (x) — a^8(x) yields a 
backscattering current 



h 
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(3) 



This contribution will be negligible at low bias V re- 
flecting the irrelevance of the p 3 -Rashba term. The im- 
portance of this nonvanishing result is, that it demon- 
strates how in principle inelastic single electron backscat- 
tering can occur at finite bias even though the perturba- 
tions Rashba SOC and electron phonon coupling preserve 
TRS. The lowest order nonvanishing matrix element for 
the oil term could be third order in ct\ which has the same 
relevance as a 3 in renormalization group (RG) sense. 
Also quadratic corrections to the linear electronic dis- 
persion which might become relevant at higher energies 
can give rise to a± backscattering. We conclude that for 
the system without Coulomb interaction the most rel- 
evant inelastic backscattering contributions allowed by 
TRS and phase space arguments cancel. This can be 
interpreted as an enhanced precision of the quantized 
conductivity of the helical edge states at finite temper- 
ature/bias going beyond the topological protection only 
pertaining to elastic scattering. 

hTLL with Coulomb interaction — Now, we want to 
investigate whether the observed robustness of the helical 
Dirac fermions against inelastic backscattering by virtue 
of the ai-Rashba SOC term, i.e. the most relevant TRS 
preserving term coupling opposite spins, persists in the 
presence of Coulomb interaction. In order to account for 
Coulomb interaction we represent the fermionic degrees 
of freedom in terms of a hTLL using the bosonization 
identity 



4'i 



(4) 



where ip± now denote the slowly varying fields with a 
factor of e kpX separated off, a is the high energy cutoff 



of the model and <p, 9 are the bosonic phase field and 
its dual respectively. In the thermodynamic limit the 
Klein factors ry* obey the algebra of Majorana fermions. 
Absorbing a prefactor of - by redefining A, the electron 
phonon coupling can be represented as 



H. 



A / dxd x (p(x) d x d(x). 



ep 



The free hTLL Hamiltonian reads 



H 



hTLL 



2 ' '''' 



n 2 Or) + 1 (dMx)) 2 



with the interaction strength parameter g < 1 modelling 
repulsive Coulomb interaction. The ai-Rashba Hamilto- 
nian with a\{x) = aS(x) in bosonized form yields 



H R = 



-?7 + r; : (d x 9(x)) cos ( \/4TTip(x) 
1 V 



where the dots denote normal ordering. 

We want to calculate the average current I(x,t) = 
-^dt(<p(x,t)) due to an applied bias V. Such expecta- 
tion values can be most easily represented for practical 
calculations in terms of the generating functional 

Z[J] = J D( v >,6,d)e lSa ~ lJ c HR+ ^ ET ' 73 ' p+ ^ jTv , (5) 

where Sq encompasses the electron phonon system with- 
out the Rashba impurity, J c is along the Keldysh contour, 
03 is a Pauli matrix in Keldysh space, and scalar prod- 
ucts like J T ip involve an integration over real space and 
time. The applied bias is modelled by E(x,t) along the 
lines of Ref. [25]. 

To make further analytical progress we now integrate 
out the phonons on the Keldysh contour. The part of the 
Lagrangian involving the phonon field d reads 



L„ 



l - ([d t df - c 2 {d x df) - \(d x <p)(d x d) 



The phonon dependent part of the action can be repre- 
sented on the Keldysh contour as 



S d = ^d T G p ^d + \d T a 3 d 2 x p. 



with the phonon propagator G p h. Performing the Gaus- 
sian integral 



Dd e iSd( - d '^ = 

in the rotated Keldysh basis 
^+\ 1 (I 1 

yields for the dissipative action 
A 2 



U(p 



(d x tf) (TlG ph (Tld x ip, 



(6) 
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/ Q QA \ 

where CTiGphCi = I nR „^ I. By this dissipative ac- 

tion the inverse free electron Green function in Fourier 
space is changed to the following dressed version 

A 2 fc 4 



{GZ\Ku>)) lptp ^{G-\k,d)) vv - 



c 2 k 2 



(7) 



This result generalizes to the Keldysh formalism a sim- 
ilar imaginary time calculation carried out in Ref. |21j . 
The retarded, advanced and Keldysh part of this Green 
function can be calculated exactly. From now on the 
free action So(ip,6) refers to the effective action where 
the phonons have been integrated out. To calculate the 
current we basically have to evaluate 



(<p(x)) 



8Z[J] 



V2 8J(r) 



(8) 



,7=0 



which can be done along the lines of Ref. [22] . The equa- 
tion of motion dt<f — —d x 9 inside the free correlators 
averaged with Sa(ip,0) remains valid for the free action 
which is dressed by the phonon dissipation. That is be- 
cause Sdiss (see Eq. (pi)) depends only on <p and not on its 
conjugate momentum d x 9. Using only this property of 
the electron phonon coupling the current to second order 
in ot\ is readily shown to be zero proving that our result 
for the case without Coulomb interaction persists in the 
nonequilibrium hTLL. 

This result can be understood on more general 
grounds. The fact that So remains quadratic in Tl v = 
d x 9 means that the following argument first brought for- 
ward in Ref. |19| for a hTLL in equilibrium without elec- 
tron phonon interaction can be used for our setup as 
well: Integrating out LI^ in the path integral representa- 
tion of the generating functional will produce terms pro- 



portional to (d t ip) , a(d t ip) cos^ 4irtp), a cos^Att ip) in 
the action. While the first contribution is the term well 
known from the free hTLL case, the second one is a pure 
gauge which can be dropped. The third term can up to 
a constant be written as cos(y/16mp) which is a two elec- 
tron TRS preserving backscattering term [TT] . Note that 
the presence of the source term J T Lp in Eq. ^ does 
not affect this argument. Thus we have found that the 
ai-Rashba term cannot lead to a single electron backscat- 
tering term in presence of any external spin independent 
dissipation which couples linearly to the electron density. 
We again point out that this result goes beyond the topo- 
logical protection of the hTLL. It is due to the quadratic 
form of So in U v , the Luttinger liquid analogue of the 
linear dispersion of helical Dirac fermions on which our 
result without Coulomb interaction relied. Bosonizing 
the a3-Rashba term implies terms up to third power in 
d x 9 thus breaking the quadratic form of So in n^. Fur- 
thermore mixed terms like (d x ip)(d x 9) 2 will occur which 
render the modifications of free Green function (see Eq. 
([T])) by the presence of phonons important. These obser- 
vations are perfectly compatible with our nonvanishing 



result for single electron backscattering in the presence 
of a 3 (see Eq.Q). 

In summary, we have studied helical Dirac fermions 
in the presence of electron phonon coupling and Rashba 
SOC, which preserve TRS. We have shown that, although 
TRS docs not provide a protection against inelastic scat- 
tering, the current carried by the helical states in pres- 
ence of a finite bias is not changed to leading order. Fur- 
thermore, we have proven that this result still holds for 
a hTLL including Coulomb interaction. The linear dis- 
persion of the helical edge states of a QSH bar has been 
nicely verified experimentally [3 and is an exact feature 
of the four band model for inverted HgTe/CdTe quantum 
wells introduced in Ref. [2J. Therefore, our analysis is 
not only interesting for the abstract model of a hTLL. It 
supports the robustness of the quantized subgap conduc- 
tance of a QSH sample beyond the well known argument 
of topological protection. Our results turn out to 
be not restricted to the coupling to longitudinal acoustic 
ID phonons. In fact, we have shown that any external 
bath coupling linearly to the electron density cannot give 
rise to inelastic single electron backscattering in presence 
of linear in k spin orbit coupling. In HgTe/CdTe layer 
structures external coupling mechanisms, e.g coupling to 
charge puddles in the bulk, are likely to cause phase de- 
coherence which gives rise to additional backscattering. 
However, such effects are not intrinsic features of the 
hTLL and can in principle be contained by improving 
the sample quality. In contrast, electron phonon cou- 
pling and Coulomb interaction are intrinsic mechanisms 
the role of which we have investigated for a generic real- 
ization of the hTLL. 
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